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Description 

[METHOD AND APPARATUS FOR 
DETERMINING SPEED AND PROPERTIES 
OF FLOWING FLUIDS USING NMR 
MEASUREMENTS^ 

Background of Invention 
[0001] Field of the Invention 

[0002] The invention relates generally to fluid characterization 

using nuclear magnetic resonance (NMR) instruments. 
[0003] Background Art 

[0004] The oil and gas industry has developed various tools ca- 
pable of determining and predicting earth formation 
properties. Among different types of tools, nuclear mag- 
netic resonance (NMR) instruments have proven to be in- 
valuable. NMR instruments can be used to determine for- 
mation properties, such as the fractional volume of pore 
space and the fractional volume of mobile fluid filling the 
pore space. A general background of NMR well logging is 



described in U.S. Patent No. 6,140,817. 
[0005] Nuclear magnetic resonance is a phenomenon occurring 
in a selected group of nuclei having magnetic nuclear mo- 
ments, i.e., non-zero spin quantum numbers. When these 
nuclei are placed in a magnetic field (B , "Zeeman field"), 

o 

they each precess around the axis of the B field with a 

o 

specific frequency, the Larmor frequency (go ), which is a 

o 

characteristic property of each nuclear species and de- 
pends on the gyromagnetic ratio 

of the nucleus and the magnetic field strength (B ) effec- 

o 

tive at the location of the nucleus, oo = vB . 

o o 

[0006] Borehole fluid sampling and testing tools such as Schlum- 
berger's Modular Dynamics Testing (MDT) Tool can pro- 
vide important information on the type and properties of 
reservoir fluids in addition to providing measurements of 
reservoir pressure. These tools may perform measure- 
ments of the fluid properties downhole, using sensor 
modules on board the tools. Alternatively, these tools can 
withdraw fluid samples from the reservoir that can be col- 
lected in bottles and brought to the surface for analysis. 
The collected samples are routinely sent to fluid proper- 
ties laboratories for analysis of physical properties that 



include, among other things, oil viscosity, gas-oil ratio, 
mass density or API gravity, molecular composition, H 2 S, 
asphaltenes, resins, and various other impurity concentra- 
tions. However, the laboratory data may not be useful or 
relevant to the reservoir fluid properties because the sam- 
ples may be contaminated by mud filtrate. 

[0007] For example, the collected fluid samples could be emul- 
sions of filtrate water and crude oil or, in wells drilled with 
oil-base muds, mixtures of reservoir crude oil and oil- 
base mud filtrate. In either case the contamination may 
render the measured laboratory data irrelevant to the ac- 
tual properties of the in situ reservoir fluids. In those cases 
where the samples brought to the surface have low or 
negligible contamination, laboratory results can still be 
tainted (e.g., by precipitation of solids caused by temper- 
ature changes). Therefore, it is desirable that formation 
fluid analysis be performed, on uncontaminated samples, 
under downhole conditions. 

[0008] Several U.S. Patents disclose methods and apparatus for 
making NMR measurements in a borehole on fluid sam- 
ples withdrawn from earth formations. U.S. Patent No. 
6,346,813 Bl issued to Kleinberg and published U.S. 
Patent Application No. 2003/0006768 by Kleinberg et al. 



disclose an NMR module on the flowline of the MDT tool 
for determining different fluid properties from magnetic 
resonance signals. The Kleinberg patent and the Kleinberg 
application are assigned to the assignee of the present in- 
vention. This patent and the application are hereby incor- 
porated by reference in their entirety. U.S. Patent No. 

6.111.408 issued to Blades etal. also discloses an NMR 
module that can be incorporated into a fluid sampling tool 
similar to the Schlumberger MDT™ tool. The NMR module 
can perform different NMR measurements including r 
(longitudinal relaxation time), T (transverse relaxation 
time), and D (diffusion constant). U.S. Patent No. 

6.111.409 issued to Edwards etal. discloses an apparatus 
and method for measuring spin-echo signals from car- 
bon-13 nuclei and a spectral method of analyzing these 
measurements to determine the ratio of aliphatic to aro- 
matic hydrocarbons from the chemical shift spectrum. 
U.S. Patent No. 3,528,000 issued to Schwede discloses 
different methods, including those using an open hole 
packer, for withdrawing fluid samples into a borehole tool 
in order to perform NMR measurements on the fluid sam- 
ples. However, none of these prior art patents disclose 
methods for determining flow speed or T of flowing flu- 



ids. 

[0009] u.S Patent 6,046,587 issued to King and Santos discloses 
methods and apparatus for measuring flow speed of mul- 
tiphase fluids flowing in a pipeline. The King and Santos 
patent teaches using the ratio of NMR free-induction de- 
cay (FID) amplitudes of signals acquired with different de- 
lay times to infer flow speed using a single NMR sensor. 
However, the FID signal is difficult to measure using per- 
manent magnets because static magnetic field variations 
will result in signals that decay too fast for reliable detec- 
tion. Moreover, the King and Santos methods for deter- 
mining flow speed do not account for the fact that there is 
a distribution of flow speeds in a pipe. In another embod- 
iment, the King and Santos apparatus consists of two sep- 
arated NMR sensors. The flow speed and fluid volumes for 
multiphase fluid flow are computed from FID measure- 
ments in the two sensors. However, the computation re- 
quires prior knowledge of the fluid T distributions. 

[0010] t wo methods for measuring the speed of fluids flowing 
through a flowline of a fluid sampling tool are described 
in a U.S. Patent Application Serial No. 10/349011, filed by 
Madio etal. on January 22, 2003. This application is as- 
signed to the assignee of the present invention and is 



hereby incorporated by reference in its entirety. Both 
methods require that the static magnetic field of the NMR 
apparatus have a gradient in the direction of the flow. The 
first method exploits the fact there is a Larmor frequency 
shift between different wait time measurements in a vari- 
able wait time pulse sequence. The measured frequency 
shifts are proportional to both the flow speed and the 
gradient. The second method is based on the fact that the 
measured phase difference between odd and even echoes 
is proportional to both the average flow speed and the 
gradient. 

[0011] while these prior art methods are useful for determining 

formation fluid properties, there remains a need for better 

methods for measuring a property of a fluid flowing in a 

formation testing tool using NMR instruments. 
Summary of Invention 

[0012] | n one aspect, embodiments of the invention relate to 

methods for determining a property of a flowing fluid by 
nuclear magnetic resonance. A method for determining a 
property of a flowing fluid by nuclear magnetic resonance 
in accordance with one embodiment of the invention in- 
cludes applying a static magnetic field to the flowing fluid; 
acquiring a suite of nuclear magnetic resonance measure- 



merits on the flowing fluid using a pulse sequence com- 
prising a spoiling pulse, a wait time, and an acquisition 
pulse sequence, wherein the suite of nuclear magnetic 
measurements have different values for the wait time; and 
fitting the suite of nuclear magnetic resonance measure- 
ments to a forward model for responses of the flowing 
fluid to derive a parameter selected from a flow speed, 
longitudinal relaxation times of the flowing fluid, and a 
combination thereof. 
[0013] Another aspect of the invention relates to methods for 
monitoring contamination in a flowing fluid being with- 
drawn into a formation fluid testing tool using nuclear 
magnetic resonance. A method for monitoring contamina- 
tion in a flowing fluid being withdrawn into a formation 
fluid testing tool using nuclear magnetic resonance in ac- 
cordance with one embodiment of the invention includes 
applying a static magnetic field to the flowing fluid; ac- 
quiring a suite of nuclear magnetic resonance measure- 
ments of the flowing fluid using a pulse sequence com- 
prising a spoiling pulse, a wait time, and an acquisition 
pulse sequence, wherein the suite of nuclear magnetic 
measurements have different values for the wait time; fit- 
ting the suite of nuclear magnetic resonance measure- 



merits to a forward model for responses of the flowing 
fluid to derive a property of the flowing fluid; and moni- 
toring a level of contamination in the flowing fluid based 
on the derived property of the flowing fluid. 

[0014] Another aspect of the invention relates to an NMR instru- 
ment for measuring a property of a flowing fluid. A nu- 
clear magnetic resonance apparatus in accordance with 
one embodiment of the invention includes a flow pipe in- 
cluding a prepolarization section and an investigation 
section, wherein the prepolarization section is upstream 
of the investigation section; a magnet disposed around 
the flow pipe for creating a static magnetic field covering 
the prepolarization section and the investigation section; 
an antenna disposed around the flow pipe at the investi- 
gation section for generating an oscillating magnetic field 
having a magnetic dipole substantially perpendicular to a 
magnetic dipole of the static magnetic field, and for re- 
ceiving a nuclear magnetic resonance signal; and a cir- 
cuitry for controlling generation of the oscillating mag- 
netic field and reception of the nuclear magnetic reso- 
nance signal by the antenna. 

[0015] other aspects of the invention will become apparent from 
the following description, the drawings, and the claims. 



Brief Description of Drawings 

[0016] FIG. 1 shows a formation fluid testing tool including an 
NMR module suitable for practicing embodiments of the 
invention. 

[0017] FIG. 2A shows a schematic of an NMR sensor in accor- 
dance with one embodiment of the invention. 

[0018] FIG. 2B shows a schematic of spin polarization profile 
along a flow pipe of an NMR sensor in accordance with 
one embodiment of the invention. 

[0019] FIG. 3A shows a diagram of a pulse sequence for acquir- 
ing variable wait time NMR measurements according to 
one embodiment of the invention. 

[0020] FIG. 3B shows a diagram of a pulse sequence including a 
spoiling pulse according to one embodiment of the inven- 
tion. 

[0021] FIGs. 4A and 4B show a schematic diagram illustrating a 
laminar flow and a non-laminar flow, respectively. 

[0022] FIG. 5 shows computed flow speeds versus input flow 
speeds in a simulation according to one embodiment of 
the invention. 

[0023] FIG. 6 shows the fitting of signal magnitudes to a suite of 
data with variable wait times for various flowing speeds 
according to one embodiment of the invention. 



[0024] FIG. 7 shows the fitting of signal magnitudes to a suite of 
data with variable wait times for various flowing speeds 
according to one embodiment of the invention. 

[0025] FIG. 8 shows computed flow speeds versus input flow 
speeds, for fluids with different T s,in a simulation ac- 
cording to one embodiment of the invention. 

[0026] FIG. 9 shows computed T versus input T for various flow 
speeds in a simulation according to one embodiment of 
the invention. 

[0027] FIG. 10 shows a flow chart of a method for determining 

the flow speed and properties of a flowing fluid according 
to one embodiment of the invention. 

[0028] FIG. 11 shows a flow chart of a method for monitoring 

contamination in a flowing fluid in a formation fluid test- 
ing tool in accordance with one embodiment of the inven- 
tion. 

Detailed Description 

[0029] Embodiments of the invention relate to methods and ap- 
paratus for determining speeds and properties of flowing 
fluids using NMR measurements. The methods disclosed 
herein for determining flow speed do not depend on static 
magnetic field gradients. 

[0030] Some embodiments of the invention also relate to meth- 



ods for measuring longitudinal relaxation time (T^of the 
fluid. The measurement of T provides a means for moni- 
toring contamination levels in the reservoir fluids that are 
being withdrawn from the reservoir into a borehole fluid 
testing/sampling tool, such as the Modular Dynamics 
Tester tool sold under the trade name of MDT™ by 
Schlumberger Technology Corporation (Houston, TX) or 
any similar fluid sampling tool. Monitoring changes in T 
with time is useful for determining when to divert or stop 
the flow in order to: (1) to make in situ NMR measurements 
of longitudinal (r ) or transverse relaxation times (T ), 
diffusion coefficients (D), or chemical shift spectra, or (2) 
collect a fluid sample that can be brought to the surface 
and sent to a laboratory for analysis. 
[0031] Some embodiments of the invention relate to methods for 
computing the viscosities of the flowing fluids. In accor- 
dance with one embodiment of the invention, the T of 

' i 

the fluid can also be used to estimate reservoir oil viscos- 
ity using known empirical correlations, provided that the 
gas-oil ratio of the fluid is known from other measure- 
ments. Another method in accordance with the invention 
may use measurements of pressure differential (AP) across 
the flowline together with the flow speed, e.g., deter- 



mined by NMR, to compute the viscosity. 
[0032] Methods of the invention may be practiced with an NMR 
instrument modified from those known in the art. The 
NMR instrument may be in the laboratory or be part of a 
wireline tool, such as a formation tester (fluid sampling 
tool) used to analyze the fluids withdrawn from the for- 
mations. 

[0033] FIG. 1 illustrates a formation fluid testing/sampling tool, 
such as Schlumberger's Modular Dynamics Testing Tool 
(MDT), which comprises several parts which enable ex- 
traction of fluids from permeable earth formations. As 
shown, the formation fluid testing tool 10 includes the 
following modules: an electronic module 11, which may 
include a processor and a memory; a hydraulic power 
module 12; a probe module 13, which may be deployed to 
make a hydraulic seal with the formation; a pumpout 
module 17; an optical fluid analyzer (OFA) 14; and a mul- 
tisample module 16. In addition, the formation fluid test- 
ing tool 10 includes an NMR sensor 15. 

[0034] A n nmr sensor 15 in accordance with an embodiment of 
the invention includes (1) a magnet (e.g., a permanent 
magnet) that is designed to produce a static magnetic 
field in the flowline (flow pipe) of the formation fluid test- 



ing tool 10, and (2) an RF antenna (which may function as 
a transmitter and a receiver) designed to radiate a oscil- 
lating magnetic field having its magnetic dipole substan- 
tially perpendicular (orthogonal) to that of the static mag- 
netic field. The frequency of the oscillating magnetic field 
is selected to be equal to the Larmor frequency of the 

1 13 

NMR sensitive nuclei (e.g., H or C) under investigation. 
Because of signal-to-noise considerations, it is preferred 
to measure 1 H nuclei in rapidly flowing fluids. For station- 
ary measurements (i.e., when fluid is not flowing), signals 
from other nuclei, including 13 C may be measured. One of 
ordinary skill in the art would appreciate that the same RF 
antenna may function as a transmitter to transmit the os- 
cillating magnetic field and as a receiver to receive the 
signals, as disclosed in U.S. Patent No. 6,346,813 Bl is- 
sued to Kleinberg. Alternatively, separate transmitter and 
receiving antennas may be used. The magnet may have a 
pre-polarization region for polarizing the fluid before it 
enters the transmitter/receiver antenna region. The pre- 
polarization may be necessary for obtaining measurable 
signals from rapidly moving nuclei. 
[0035] FIG. 2A shows an illustration of an NMR sensor 15 in ac- 
cordance with one embodiment of the invention for deter- 



mining flow speed and other properties of a fluid flowing 
in a flowline 22. As shown, the sensor 15 includes a non- 
conductive and non-magnetic flowline 22, a permanent 
magnet 24, and an antenna 26 for transmitting and re- 
ceiving radio frequency (RF) signals. 
[0036] As shown, the flowline (flow pipe) 22 includes a prepolar- 
ization section 28 that is upstream of an investigation 
section 29. The magnet 24 is disposed around both the 
prepolarization section 28 and the investigation section 
29. The magnet 24 may be a permanent magnet or an 
electromagnet magnet. As shown, the flowline 22 is made 
of a non-conductive and non-magnetic material, such as 
a composite or polymer material. If the flowline 22 is 
made of a conductive or magnetic material (e.g., steel), 
then the antenna 26 is preferably located inside the flow- 
line. 

[0037] For the sensor shown in FIG. 2A, the antenna (solenoid 

coil) 26 functions as both a transmitter and a receiver. Al- 
ternatively, the NMR sensor 15 may include two antennas: 
one functions as a transmitter and the other as a receiver. 
Furthermore, the antenna may comprise a solenoid coil, a 
loop, or a saddle coil. The antenna 26 substantially covers 
the investigation section 29 of the flowline 22. 



[0038] The NMR measurements according to methods of the in- 
vention includes a suite of variable wait time (W) measure- 
ments. Prior to each wait time (W), the magnetization is 
first spoiled by pulses designed to "kill" the magnetization 
so that M =M =M = 0. Following the spoiling pulse 

x y z 

and the wait time, a 90-degree pulse followed by a 
180-degree pulse (e.g., a spin-echo pulse) is applied to 
the transmitter to generate a spin echo. 

[0039] The receiving antenna measures the in-phase and 

quadrature components of the spin-echo signal. The 
measurement is repeated for a number of different wait 
times. In other embodiments of this invention, multiple 
180-degree pulses may be applied to produce multiple 
spin-echo signals. In addition, one of ordinary skill in the 
art would appreciate that spin-echoes may also be ob- 
tained using pulses other than 90- and 180-degree 
pulses, for example, using composite pulses. 

[0040] The amplitudes of the detected spin-echo signals for the 
different wait times depend on flow speed, wait time, re- 
ceiver and transmitter lengths, magnet pre-polarization 
length, and the T distribution of the fluid. All of these 
parameters, except the flow speed and the r distribution 
of the fluid, are either fixed by the sensor design or by the 



pulse sequence. If sufficient sets of measurements are 
available, these parameters may be derived by fitting the 
measured signals to a proper model that simulates the 
NMR response of the flowing fluid. That is, the data may 
be interpreted by forward modeling. 
[0041] The method disclosed below employs a theoretical for- 
ward model that allows the computation of the detected 
echo signals for any combination of the aforementioned 
pulse sequence, sensor parameters, flow speed, and r 
distribution. The forward model may be fit, by inversion, 
to the measured echo signals to determine the apparent 
flow speed and T distribution. Alternatively, if the flow 
speed is known from other measurements then the for- 
ward model can be fit by inversion to determine the T 

7 1 

distribution. 

[0042] Sensor and Variable Wait Time Pulse Sequence 

[0043] FIG. 2B shows that the speed and position dependent po- 
larization profile (/(v,z)) over the length (/ ) of the receiver 

a 

coil (shown as 26 in FIG. 2A) is made up of two parts - 
one part comes from fresh spins that are "pre-polarized" 
as they travel in the static magnetic field into the receiver 
during the wait time (W) and the second part that comes 
from spins that were in the receiver coil and are "re- 



polarized" during the wait time (W). The pre-polarization 
length (/ ) shown in FIGs. 2A and 2B corresponds to the 
length of the prepolarization region 28. 
[0044] The position and speed dependent polarization function is 
given by the equation, 

f(y,x) =l-exp(-^^. CO 

where 7tz,v) is the polarization time for a spin with posi- 
tionzand speed v, and T is the longitudinal relaxation 
time of the spins. The polarization time, 7tz,v), is defined 
as: 

L + z 

T(z,v)=-Z for 0<x<vW (2$ 

V 

and, 

Y{s,v)=W for v W <2<l u . (2b) 

[0045] it is clear from the above equations that the speed (v) and 
position (z) dependent polarization function, f(v,z), also 
depends on / , w, and T ^ However, these dependencies 
are implicit in Eq. 1 because polarization time, T(z,v), is 
used to simplify the notation. These equations can be 
more easily understood following a discussion of the vari- 
able wait time (VWT) pulse sequence shown in FIGs. 3A 
and 3B. 



[0046] FIG. 3A shows a suite of variable wait time pulses consist- 
ing of N measurements. N is preferably on the order of 
10. After each measurement in the suite, one or more 
spoiling pulses (collectively denoted by the pulse S) are 
applied to destroy any remnant magnetization before 
starting the wait time for the next measurement. The du- 
ration and frequencies of the spoiling pulses are instru- 
ment dependent and can be determined empirically. 

[0047] As s hown schematically by the "tadpole" in FIG. 3B, fol- 
lowing each wait time a 90-degree excitation pulse is ap- 
plied to rotate the longitudinal magnetization onto the 
transverse plane. The signal from the transverse magneti- 
zation rapidly dephases due to inhomogeneity (or other 
factors) in the static magnetic field, but is refocused by 
the 180-degree pulse to produce a spin-echo. After the 
spin-echo signals are acquired, a spoiling pulse is applied 
to remove the magnetization of the spins within the re- 
ceiver coil (shown as 26 in FIG. 2A). Those skilled in the 
art will recognize that a variable wait time sequence like 
the one shown in Fig. 3 can also be used observe the 
free-induction decay signals instead of a spin echo. That 
is, the 180-degree pulse shown in Fig. 3B can be omitted 
and the free-induction decay signal can be observed fol- 



lowing the 90-degree pulse. Only minor and obvious 
modifications to the forward model are required. Those 
skilled in the art will also recognize that other T - 
sensitive sequences, such as inversion recovery or satura- 
tion recovery pulse sequences, could also be employed 
with obvious modifications to the forward model. Both the 
inversion-recovery and the saturation-recovery pulse se- 
quences are generally referred to as ""^-relaxation inves- 
tigation pulse sequences" in this description. 
[0048] The pulse sequence used to acquire the NMR signals (FID 
or echoes) are generally referred to as an "acquisition 
pulse sequence" in this description. An acquisition pulse 
sequence may include a single 90-degree pulse, a spin- 
echo pulse (i.e., a 90-degree pulse followed by a 
180-degree pulse), and the variants of the spin-echo 
pulse that include multiple 180-degree (refocusing) 
pulses following the 90-degree pulse (e.g., a CPMG). 
Thus, the "variable wait time pulse sequence" as shown in 
FIGs. 3A and 3B comprises a spoiling pulse, a wait time, 
and an acquisition pulse sequence. The spin-echo pulse 
sequence and its variant (e.g., CPMG) are generally re- 
ferred to as a "spin-echo pulse sequence" in this descrip- 
tion. That is, a "spin-echo pulse sequence" not only com- 



prises a single 90-degree pulse and a single 180-degree 
(refocusing) pulse, but also comprises multiple 
180-degree (refocusing) pulses after the 90-degree pulse. 
[0049] During each wait time, fresh or pre-polarized spins move 
into the antenna region. Eq. 2a shows that in the region, 

0 < z < vW, 

fresh spins have entered the antenna during the wait time. 
The length of this region depends on both the flow speed 
(v) and the wait time (W), i.e., it depends on the product of 
the flow speed (v) and the wait time (W). The polarization 
time for these fresh spins is independent of the wait time. 
Instead, it depends on the duration that the spins have 
been exposed to the static magnetic field since they en- 
tered the permanent magnet. This is because this portion 
of the fluid was outside the transmitter/receiver antenna 
(shown as 26 in FIG. 2A) when the spoiling pulse is ap- 
plied and the magnetizations of the spins in this portion 
of the fluid were not "killed" by the spoiling pulse. On the 
other hand, as can be seen from Eq. 2b, the re- 
polarization of the spins in the adjoining region of length, 
/ - vW, is controlled by the wait time (W) because these 

a 

spins were in the transmitter/receiver region when the 
spoiling pulse was applied. The magnetizations of these 



spins were removed by the spoiling pulse, and any polar- 
ization detected by the receiver is due to repolarization 
during the wait time (W). If the wait time is long or the 
flowing speed is fast enough such that the product, vW, 
exceeds the antenna length, then only fresh spins (those 
that enter the receiver antenna after the spoiling pulse) 
are measured and the polarization function is indepen- 
dent of w. 

[0050] Derivation of the NMR Forward Model for Flowing Fluids 

[0051] The methods disclosed herein preferably employ forward 
models for predicting the responses of the NMR sensor 
for a fluid with a given flow speed and distribution of T 
relaxation times. The response of the sensor also depends 
on the wait times and physical parameters of the sensor 
such as the antenna length, the magnet pre-polarization 
length, and the radius of the flowline. For a given wait 
time and sensor design the only variables in the forward 
model are the flow speed and the T distribution of the 
fluid. The flow speed and r distribution are preferably 
determined by inversion of a forward model. The forward 
model is derived in the following paragraphs. 

[0052] jo accurately model the NMR responses of flowing fluids, 
the fact that the velocity profile of a laminar flow for a vis- 



cous fluid flowing in a pipe is parabolic should be taken 
into account. See e.g., Victor L. Streeter, "Fluid Mechanics" 

th 

McGraw-Hill Book Co., 5 Edition, p. 244. In a laminar 
flow, the maximum flow velocity, v , occurs at the axis of 

m 

the pipe, while the velocity is zero at the wall of the pipe. 
The laminar flow regime in circular pipes is characterized 
by the values of Reynolds number, 

R 12000 -3000, 

where the exact upper limit for laminar flow depends on 
the surface roughness of the pipe. The Reynolds number 
is defined by, 



where r is the radius of the flowline (pipe), v is the aver- 

o 

age flow speed, □ is the fluid mass density, and nis the 
viscosity of the fluid. In contrast to a laminar flow, turbu- 
lent flow has a chaotic component and is much more dif- 
ficult to model. One feature of turbulent flow is a flatten- 
ing of the velocity parabolic profile. FIG. 4 shows the ve- 
locity profiles for laminar flow (A) and non-laminar flow 
(B) in a circular pipe. To develop a forward model, an ex- 
pression for the echo amplitudes measured by the vari- 
able wait time pulse sequence shown in FIGs. 3A and 3B is 



needed. If A represents the amplitude of a phase- 

i 

corrected echo signal measured with a wait timePf ., then 
A . can be written in the form of a volume integral over the 
sensitive volume (Q) of the receiver antenna, i.e., 

[0053] The inner integral in Eq. 4 is a T ^weighted integral over 
a speed and position dependent polarization function. The 
weighting function is the T distribution of the fluid, P(T 
). The volume integral is weighted by the spatial sensitiv- 
ity function, S (r), that can be determined by measure- 

a 

ments or calculations. To simplify the following discus- 
sion, it is useful to consider a receiver having a uniform 
spatial sensitivity so that, 

where r is the radius of the flow pipe. Using Eq. 5 and 

o 

cylindrical coordinates (r,n, z)the integral in Eq. 4 can be 
written in the following form: 



[0054] | n arriving at Eq. 6, the angular integral was performed 
assuming the sensitivity function has azimuthal symme- 
try. Note that the polarization function in the integrand of 



Eq. 6 is given by Eqs. 1 and 2a- 2b. To proceed with the 
evaluation of A ,the speed dependent polarization func- 

i 

tion may be obtained by integrating over the z co- 
ordinate, i.e., 

The above integrals can be performed analytically to pro- 
duce: 

(8 a) 

for w -v < / , and for w -v>l , the integral is given by: 

i a i a 

/(vM,Ii,/ J 8^=1 -JLL-a-eKpC-A-JJ-eKpC-i.) . (8b) 

[0055] Before proceeding further, it is useful to make some ob- 
servations in order to provide some insights into Eqs. 
8a-8b. The first two terms in Eq. 8a account for the pre- 
polarization of the spins, whereas the last term is due to 
re-polarization of spins that are not removed by outflow 
from the receiver antenna region during the wait time (W). 
Note that Eq. 8b can be obtained from 8a by setting W . -v 
= / . Once the condition w -v>l is met, the sensitive re- 

a i a 

gion is totally filled with fresh spins and the polarization 
function is independent of w .. 



[0056] The functions in Eqs. 8a and 8b represent the radial sen- 
sitivity functions after integration over the z co-ordinate. 
Next, it is necessary to integrate these functions over the 
radial co-ordinate. The radial dependence of the radial 
sensitivity function arises from the radial variations of the 
flow speed in the flowline (see FIG. 4). For laminar flow 
the velocity profile is parabolic and can be written in the 
form, 



where r is the radius of the flowline and v is the maxi- 

o m 

mum flow speed on the axis of the flowline (i.e., at r = 0) 
as shown in FIG. 4. It follows from Eq. 9 that the average 
flow speed with a laminar flow is v 12. While Eq. 9 de- 



scribes a commonly used model for laminar flow, an alter- 
native laminar flow model may be described as follows: 



where n is typically between 5 and 10. Embodiments of 
the invention can use either of these two laminar flow 
models or any variation thereof. For clarity, the following 
description is based on the model described in Eq. 9. 




m 




However, one of ordinary skill in the art would appreciate 
that the model described in Eq. 9a can substitute for the 
model of Eq. 9 in the following description. 
[0057] Consider the following integral over the radial co- 
ordinate, 

F^Jul^lM =\ll"rf{<rXTi,l p ,l ± ,Wi)dr (10) 

[0058] Because Eq. 9 describes the relationship between the ra- 
dius rand the flow speeds v, the integral over the radiusr 
shown in Eq. 10 can be transformed into an integral over 
flow speeds vby using Eq. 9 to make a change of vari- 
ables, e.g., 

dv=-^.v m dr. (11) 

Eqs. 9 and 11 may be used to change the integration vari- 
able and limits in Eq. 10 to give, 

It is useful to make another change of variables in the in- 
tegral in Eq. 12. Making the change of variables 

in the integral in Eq. 12 leads to the integral, 



(13) 



which is in a form that can be accurately performed nu- 
merically using the Gaussian integration method. Assem- 
bling the pieces, Eq. 6 for the echo amplitudes can be 
rewritten in the following form: 

A=l^Ji)^WM^i- (14) 
Note that with known values of/ , / , and W , the model 

p a i 

echo amplitudes depend only on P(T ), the T distribu- 
tion, and the maximum flow speed, v . As described in a 

m 

later section, estimates of P(T ) and v may be obtained 

1 m 

by using suites of variable wait time data to invert the 
model. For laminar flow, v and the flowline radius r 

m o 

completely specify the velocity profile in Eq. 9. 

[0059] The expression in Eq. 13 was derived based on the 

parabolic velocity fluid profile described in Eq. 9. This is 
appropriate for laminar flow. However, for non-laminar 
profile shown in FIG. 4, the equation derived in Eq. 13 is 
not valid and should be modified. 

[0060] The velocity profile for the non-laminar flow depicted in 
FIG. 4 can be described by the following equations: 

vtr) = v n for 0<r<a, (13a) 



and, 

Using Eqs. 15a and 15b, the integral in Eq. 10 can be ex- 
pressed as the sum of two integrals, i.e., 

The first integral on the right is straightforward. Trans- 
forming the integration over r to integration over v by us- 
ing Eq. 15b performs the second integral. The last step is 
to introduce the change of variables 

to convert Eq. 16 into the following form: 

Fiv m ^,T l J p ,L,W ! )=^ r f(v m J 1 J p J A ,W l ) +^-^-J p 1 /(v^IUp,/^)^. (17) 

Note that if a= 0, Eq. 17 reduces to Eq. 13 for laminar 
flow. Moreover, the model echo amplitudes in Eq. 14 now 
depend on the parameter a in addition to their depen- 
dence on v and P(T ). By inversion of the forward 

m 1 

model, it is feasible to estimate the parameter, a, and, 
therefore, to infer the onset or existence of turbulent or 
non-laminar flow. 



[0061] Numerical Integration 

[0062] The integral in Eq. 17 can be expressed as a weighted 
sum of the integrand over a set of fixed abscissas, i.e., 

Fiv^aJ^lyJM =^f(v m J 1 J p J a ,W l ) + b*^^£ ^f^^hJ^M. (IS) 

The weights, w , and the abscissas, x ,for nth-order 

k k 

Gaussian integration can be found in tables given in 
"Handbook of Mathematical Functions With Formulas, 
Graphs, and Mathematical Tables," edited by M. 
Abramowitz and L. A. Stegun, National Bureau of Stan- 
dards Applied Mathematics Series 55, June 1964. For ex- 
ample, in the table below the abscissas and weights for 
6th order Gaussian integration are shown. 



Abscissas and Weights for 6th Order G aussi an Integration 






0.233619186083197 


0.467913934572691 


0.66120938646(5265 


0.360761573048139 


0.932469514203152 


0.171324492379170 



The abscissas given in the table are for integrals in the 
range from (-1,1). The abscissas in Eq. 18 can be ob- 
tained from the latter using the equation, 



(19) 



In performing the summation in Eq. 18, it is useful to re- 
call that the "integrand" is the speed dependent polariza- 
tion function in Eqs. 8a or 8b depending on whether 

lafvnWi or In A? m fflj . 

[0063] inversion of the Forward Model 

[0064] Once the NMR data of the flowing fluid are acquired with 
variable wait times, they are fitted to the forward model to 
provide the flow speed and T information of the flowing 
fluid in accordance with embodiments of the invention. 
The inversion involves varying the parameters v , P(T ), 

m 1 

and a in the forward model until the measured echo am- 
plitudes fit the ones predicted by the model. There are 
many different ways of fitting the measured data to the 
model. The following describes two approaches as exam- 
ples. One of ordinary skill in the art would appreciate that 
other approaches may also be used without departing 
from the scope of the invention. 
[0065] inversion Method No. 1 

[0066] Measured NMR data consist of signals measured in the in- 
phase and quadrature channels of an NMR instrument. 



The quadrature signals are ninety degrees out of phase 
with respect to the corresponding signals detected in the 
in-phase channel. 
[0067] The measured signals include noises. Thus, to fit the 
model to the measurement data, some measurement 
noises should be added to the signals derived from the 
model. The following equations illustrate this approach: 

(20b) 

[0068] W here 

K and A l A 

are measured amplitudes in the in-phase and the quadra- 
ture channels, respectively, for the measurement with wait 
time, w ; A is the model predicted signal amplitude; □ is 

i i 

the model predicted signal phase; and I and I ,are 

i,r i,q 

noises in the signal and quadrature channels, respectively. 
I and I , are assumed to be zero-mean Gaussian 
noises with equal variances, a in the two channels. Using 
Eq. 14, the signal amplitude A . can be written in the fol- 
lowing form: 



J-l 



(21) 



where the integral over r in Eq. 14 has been replaced by 
a summation over theiV c amplitudes (a ( ) that are used to 
approximate the T distribution. The following quantities 
are defined: 

F^v^ai^Fiy^Y^XJtW). (22) 

and recalled from Eq. 18 that includes the parameter, a, to 
describe the radius over which the parabolic front is flat- 
tened out because of the non-laminar flow. The discrete 

set of relaxation times, T , are fixed and can be se- 

i, i 

lected to be equally spaced on a logarithmic scale. Be- 
cause the in-phase and the quadrature signals are 90 de- 
grees out of phase to each other, the signal phase, may 
be estimated from the measured in-phase and quadrature 
amplitudes using the following equation: 



Lie 

= ™*ai<£4 ), (23) 

24, 



where the caret over □ denotes the fact that it is an esti- 
mate; and N is the number of measurements (different 

m 

wait times) in the variable wait time measurement suite. 
The phase estimate obtained from Eq. 23 may then be 



used to compute the phase-corrected signal plus noise 
amplitudes and the root mean square (rms) noise accord- 
ing to the method disclosed in U.S. Patent No. 5,291,137 
issued to Freedman ("the Freedman patent"). This patent 
is incorporated by reference in its entirety. Briefly, 

= \ r cos $ + 4 <if sin 4> (24a) 

and, 

3f-> =4 r sW-4*<W. (24b) 
The Freedman patent shows that, 

<^ + J> = 4 (25) 

where the angular bracket denotes the expectation value 
andA is the model predicted amplitude as described in Eq. 

i 

14. The rms noise estimate can be computed from the 
equation, 



The measured amplitude data may be fit to the model 
amplitudes by minimizing the sum of the squared fit er- 
rors, for example, by minimizing an objective function: 



with respect to the parametersv , a, and the set of ampli- 

m 



tudes a . The last term in Eq. 27 is used to stabilize the 
solution by reducing the effects of noise amplification that 
can occur in under determined inverse problems. The pa- 
rameter y is a regularization parameter and y> 0. 

[0069] N 0 te that "spin-dynamics" effects in the magnetic field 
gradient of a permanent magnet can lead to an attenua- 
tion of the measured echo amplitudes. See Hurlimann and 
Griffin, "Spin-Dynamics of Carr-PurceU-Meiboom-GiU-like Se- 
quences in Grossly Inhomogeneous B and B Fields and Applica- 
tion to NMR Well Logging" J. Magn. Reson. v. 143, 120-135 
(2000). Therefore, the echo amplitudes in Eq. 27 are cor- 
rected for spin-dynamics attenuation. 

[0070] The minimization in Eq. 27 can be performed using a 

number of commercially available optimization programs 
that can minimize linear or non-linear functions subject to 
constraints. One such program is the double precision 
subroutine called "DLCONG" that is available as part of a 
library of optimization codes sold under the trade name of 
IMSL™ by the Visual Numerics Corp. (San Ramon, CA). The 
DLCONG code is more robust if analytical functions can be 
provided for the partial derivatives of the objective func- 
tion. The derivatives for the objective function can be 
computed analytically. For example, from Eq. 27, the fol- 



lowing derivatives may be obtained, 

™=-2-Z(.A i W-^,a,{<i l }))F^ Vm , a ) +2-ya k , (28) 



and, 



J-1 *J 



and, 



The derivative in Eq. 29 of the objective function requires 
differentiation of P. with respect to v According to Eq. 

i,l m . 

17, this differentiation requires the derivative of the 
speed-dependent polarization function, i.e., 

a 2 IfutvJjT.-a) 1 ^ 3/,j(v m jj 



In taking the derivatives in Eq. 31, the order of integration 
and differentiation is interchanged and a compact nota- 
tion for the speed-dependent polarization functions is in- 
troduced as follows: 

f.Av n -®=fM,Tv,I.,I„Wn, (32) 



which can be obtained from Eqs. 8a and 8b by making the 
substitution, 



. Because the inversion depends critically on the accuracy 
of these derivatives, it is useful to display the detailed re- 
sults. The derivatives in the summation of Eq. 31 are 
given by, 

Jv m L v m l* v m s-h.i T i.i 

C33) 

when 
and 

(34) when 

Note that the partial derivatives for the first term on the 
right hand side of Eq. 31 can be obtained from Eqs. 33 
and 34 by setting 1. The partial derivative of F. ( with 
respect to a in Eq. 30 is straightforward and can readily be 
obtained by differentiation of Eq. 17. 
[0071] inversion Method No. 2 

[0072] The approach described in the previous section fits the 

phase corrected data to the relaxation model. An alterna- 



tive approach described below simultaneously fits the two 
channel data to the forward model. An error or objective 
function, similar to that of Eq. 27, may be defined as: 

l-l i-l M 

(35) 

where 

4 lt =4sW , (36b) 

and 

= ^ sin ^ . 

The objective function of Eq. 35 may be minimized with 
respect to v , a, and the two sets of amplitudes {a } and 

m l,r 

{a }. The amplitudes in the T distribution may be com- 

l,q 1 

puted directly from the two sets of amplitudes using the 
equation, 



The signal phase can also be estimated from these ampli- 
tudes using the equation, 



2^ 



i-1 



(38) 



where AMs the number of components used to approxi- 
mate the T distribution. Again, the minimization of the 
objective function in Eq. 35 may be accomplished with any 
of a number of commercially available programs (e.g., DL- 
CONG) and the partial derivatives discussed in the previ- 
ous section. 
[0073] Numerical Simulations 

[0074] Embodiments of the invention can accurately estimate the 
flow speed and T values (either a distribution or an aver- 
age of the T values of the fluid) in a flowing fluid. The 
following describes several computer simulation experi- 
ments that demonstrate the accuracy with which v and T 

m 

can be estimated using embodiments of the invention. 
The simulations were performed for an NMR sensor (as 
shown in FIG. 2) with / = 0.5 in., / = 10.0 in. and r = 

a p o 

0.1 in. There were 35 wait times in the variable wait time 
suite ranging from 0.001 s to 0.72 s. The average flow 
speeds used in the simulations were 2.5, 5.0, 10.0, 15.0, 
20.0, 25.0, 30.0, 35.0 and 40.0 in/s. 
[0075] FIG. 5 shows a comparison of the average fluid velocity (= 



v /2) versus the estimated average fluid velocity. The 

m 

simulation was performed assuming an rms noise ampli- 
tude of 0.1 p.u. This is a relatively high signal-to-noise 
ratio considering that the signal amplitude is 100 p.u. 
when measured on a fully polarized stationary fluid. The 
simulation shown in FIG. 5 was performed for three differ- 
ent fluids having a range of T values. The T distribu- 
tions used in the simulations were single valued (i.e., 
Dirac delta functions). Note that there is good agreement 
for all three fluids between the estimated average flow 
speed (y-axis) and the input values (x-axis) used in the 
simulation. There is some loss of accuracy at the higher 
average flow speeds (e.g., above about 30 in/s). However, 
the estimates are reasonably good considering the high 
flow speeds. 

[0076] FIG. 6 shows the variable wait time data suites as a func- 
tion of flow speed, for a fluid having r = 0.316 s, and 
the post-inversion fits of the forward model to these data. 
Note that in FIG. 6 the measured echo amplitudes de- 
crease significantly with increasing flow speed due to in- 
complete polarization. This may account for the loss of 
accuracy in the flow speed estimates in FIG. 5 at higher 
flow speeds. This effect is exacerbated for fluids with long 



[0077] FIG. 7 shows the variable wait time data and the post- 
inversion fit for a fluid having T = 3 s at various flow 
speeds. The quality of the fit is assessed by a dimension- 
less chi-squared "goodness of fit" parameter defined as, 

JV W a- m 

2 

For a perfect fit, x = 1, to within noise errors. 

[0078] N 0 te the significant reduction in the echo amplitudes with 
increasing flow speed. This translates into a significant 
reduction in the measurement signal-to-noise ratios. An- 
other important feature of the data in FIG. 7 is the com- 
pression of the data (i.e., the curves are closer together) 
as the flow speed increases. Thus, the relative sensitivity 
to flow speed of measurements made at high flow speed 
is reduced. Both of these effects, reduced signal-to-noise 
ratio and reduced sensitivity, mean that it is more difficult 
to accurately invert a variable wait time data suite at 
higher flow speeds. 

[0079] FIG. 8 shows the results of the same simulation as in FIG. 
5, but with an rms noise of 1.0 p.u. instead of 0.1 p.u. A 
comparison of FIGs. 5 and 8 reveals that the quality of the 
flow speed estimation is significantly degraded by a ten- 
fold reduction in the signal-to-noise ratio. This clearly 



shows the importance of having high signal-to-noise ratio 
measurements. 
[0080] For fluids with broad r distributions, the geometric 

mean of the distributions or other moments of the distri- 
bution can be used to monitor and characterize the fluids 
in the flowline. 

[0081] Estimation of Fluid Viscosity from the Measurements 

[0082] a s noted above, some embodiments of the invention re- 
late to methods for estimating fluid viscosity from NMR 
measurements of a flowing fluid. Once the flow speed of 
the fluid is determined from a suite of NMR measurements 
having variable wait times according to methods de- 
scribed above, the viscosity of the fluid can also be esti- 
mated according to methods described in the following 
sections. Two methods will be described for estimating 
fluid viscosity (□) using the variable wait time measure- 
ments made on flowing fluids in the flowline. One of ordi- 
nary skill in the art would appreciate that these two ap- 
proaches are for illustration only, and other methods may 
be used without departing from the scope of the inven- 
tion. 

[0083] The first method computes an estimate of the viscosity (□) 
using the average flow speed (v) estimated from the NMR 



measurements and the measured pressure drop, □ P, 
across a length L of the flowline. The pressure drop, □ P, 
may be obtained from pressure sensors equipped on the 
NMR sensor or the formation testing tool. The computa- 
tion of fluid viscosity from the flow speed (v) and the 
pressure drop (□ P) may be accomplished according to the 
well-known fluid mechanics formula (See e.g., Victor L. 

th 

Streeter et al., Fluid Mechanics, 5 Edition, p. 244, McGraw- 
Hill, NY.): 

where the average speed v of the fluid is equal to v 12. 

m 

Because L and r are defined by the apparatus, Eq. 40 may 

o 

be converted to a more general form, 

17-*- — , (41) 

v 

where the proportionality constant K may be determined 
empirically from measurements made on fluids with 
known viscosities measured at different speeds with 
known pressure gradients. Using v, □ and an independent 
measure of the mass density of the fluid, it is possible to 
compute the Reynolds number, R, from Eq. 3. The compu- 
tation in Eq. 40 or 41 is valid for any viscous fluid and is 



not restricted to crude oil or brine. 
[0084] The second approach for estimating crude oil viscosity is 
to use the geometric or logarithmic mean of the T distri- 
bution. The logarithmic mean of the T distribution is re- 
lated to the oil viscosity by a known correlation (See e.g., 
Freed man et al, "A New NMR Method of Fluid Characterization in 
Reservoir Rocks: Experimental Confirmation and Simulation 
Results"SPE Journal, December 2001, 452-464), 

*' = ^W (42) 

where □ is the oil viscosity in centipoise, k is an empiri- 

o 

cally determined constant for crude oils, ris the tempera- 
ture in Kelvin, T is the logarithmic mean of the oil T 

1, LM 1 

distribution and /(GOR) is an empirically determined func- 
tion of the solution gas-oil ratio (COR). The COR in Eq. 42 
may be determined from other measurements in the flow- 
line or might be known from local or prior knowledge of 
the reservoir fluid properties. 
[0085] The computed T distribution or the viscosity of the flow- 
ing fluid computed from Eq. 40, 41, or 42 can be used to 
monitor changes in the contamination of the fluid in the 
flowline. This application will be described in the following 
sections. 

[0086] FIG. 10 illustrates a method 100 in accordance with em- 



bodiments of the invention. In accordance with this 
method, a suite of NMR data are acquired with variable 
wait times (step 102). The pulse sequences for the data 
acquisition may include a spoiling pulse, a wait time, and 
an acquisition pulse sequence, as shown in FIGs. 3A and 
3B. The acquisition pulse sequence may include a simple 
90-degree excitation pulse, a spin-echo pulse 
(90-delay-180), a Carl-Purcell-Meiboom-Gill (CPMG) 
pulse, or any variant of the spin-echo or CPMG pulse. A 
single spin-echo pulse is shown in FIG. 3B for illustration. 
The 90-degree and the 180-degree pulses shown may be 
replaced with other pulses (e.g., composite pulses) with- 
out departing from the scope of the invention. 
[0087] a s noted above, in addition to the variable wait time pulse 
sequence shown in FIGs. 3A and 3B, embodiments of the 
invention may also use other T investigation pulse se- 
quences. For example, an inversion-recovery pulse se- 
quence or a saturation-recovery pulse sequence may be 
used instead. If a saturation-recovery pulse sequence is 
used, the model fitting process is similar to that using the 
spoiling pulse sequence. However, if an inversion-recov- 
ery pulse sequence is used, then the spin polarization 
function of Eq. 1 would be modified. 



[0088] Once the suite of variable wait time NMR measurements 
are acquired, they may be fit to a forward model using ei- 
ther of the inversion methods discussed above (step 104). 
The inversion process produces the flow speed (v) and the 
T values of the flowing fluid. The T values thus obtained 
may be a T distribution or a logarithmic mean of the T 
distribution. 

[0089] The flow speed (v) and/or the T values obtained from fit- 
ting the NMR data to a flow model may be used to com- 
pute the viscosity of the flowing fluid (step 106). As noted 
above, the estimated flow speed (v) together with an inde- 
pendently measured pressure drop (AP) across the length 
L of the NMR sensor may be used to compute the fluid 
viscosity (n) according to Eq. 40 or 41. Similarly, the esti- 
mated T values or its logarithmic mean may be used to- 
gether with an independently determined COR function, 
f(GOR), to compute the fluid viscosity (n), as shown in Eq. 
42. However, any variations of these approaches may be 
used. 

[0090] As noted above, embodiments of the invention may be 
used to monitor fluid contamination as they are with- 
drawn from the formation. With the ability to compute the 
T and viscosity of the flowing fluid according to methods 



shown in FIG. 10, it becomes possible to monitor the pu- 
rity of the formation fluids as they are drawn into the for- 
mation fluid testing tool. 

[0091] FIG. 11 shows a method 110 for such monitoring in ac- 
cordance with a method of the invention. According to 
this method, formation fluids are withdrawn into the test- 
ing tool (step 112) and the viscosity or T values of the 
fluids are monitored (step 114). Because the early portion 
of the withdrawn fluids are likely contaminated with mud 
filtrates, their properties (e.g., viscosity and T values) are 
likely different from un-contaminated formation fluids. As 
the flow of the formation fluids continues, the level of the 
mud filtrate contamination will gradually decrease. There- 
fore, by monitoring the viscosity or the T values of the 
flowing fluids, it is possible to discern when the mud fil- 
trate contamination diminishes to an insignificant level, 
for example by monitoring the viscosity or T values ap- 
proaching a steady state value that reflects the properties 
of the pure formation fluids. 

[0092] By querying whether the viscosity or the T values reach a 
stable value (step 116), a decision can be made whether 
to continue the withdrawal of the formation fluids. For ex- 
ample, if the viscosity or T values has not stabilized, the 



withdrawal process and the monitoring will be continued. 
If the viscosity or the T values stabilize, then the with- 
drawn fluid samples may be saved for later analysis, or a 
downhole measurement may be performed in situ (step 
118). 

[0093] Various NMR measurements have been performed down- 
hole. For example, U.S. Patent No. 6,229,308 issued to 
Freedman and U.S. Patent No. 6,462,541 Bl issued to 
Venkataramanan et al. disclose NMR measurements under 
downhole conditions. These measurements can provide 1^ 
(transverse relaxation times) distribution, D (diffusion) 
distributions of the formation fluids, T -D maps, fluid 
typing, and fluid viscosity (n). 

[0094] Method 110 shown in FIG. 11 is functionally equivalent to 
many monitoring techniques currently used to monitor 
contamination in fluids withdrawn by a formation fluid 
testing tool. Other methods may monitor, for example, 
color, optical density, or GOR of the fluids. Thus, method 
110 shown in FIG. 11 provides an alternative to these ex- 
isting methods. 

[0095] The methods as shown in FIGs. 10 and 11 may be imple- 
mented as computer programs and stored on a computer 
or a processor used with well logging, for example the 



processor or circuitry included in the electronic module 11 
shown in FIG. 1. 
[0096] Advantages of the invention include the following. Em- 
bodiments of the invention provide convenient methods 
for estimating flow speed, T values, and viscosity of a 
flowing fluid. Methods of the invention may be used on a 
formation fluid testing tool or in any other type of appa- 
ratus and application where one needs to monitor fluid 
flow. If used on a formation testing tool, methods of the 
invention may be used to monitor contamination of the 
formation fluids as they are removed from the formation. 
These methods ensure that formation fluids substantially 
free from contamination are removed for analysis. By pro- 
viding methods for monitoring contamination in formation 
fluids while they are being withdrawn, methods of the in- 
vention can alleviate problems associated with contamina- 
tion in formation fluid analysis. In addition, methods of 
the invention also provide NMR measurements downhole, 
after the monitoring indicates that contamination is in- 
significant. Being able to measure the fluid properties un- 
der conditions close to reservoir conditions can provide 
fluid properties that are representative of the in situ fluid 
properties. 



[0097] while the invention has been described with respect to a 
limited number of embodiments, those skilled in the art, 
having benefit of this disclosure, will appreciate that other 
embodiments can be devised which do not depart from 
the scope of the invention as disclosed herein. Accord- 
ingly, the scope of the invention should be limited only by 
the attached claims. 



